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We calculate the power spectrum of metric fluctuations in inflationary cosmology starting with 
initial conditions which are imposed mode by mode when the wavelength equals some critical length 
£ c corresponding to a new energy scale M c at which trans-Planckian physics becomes important. 
In this case, the power spectrum can differ from what is calculated in the usual framework (which 
amounts to choosing the adiabatic vacuum state). The fractional difference in the results depends 
on the ratio ao between the Hubble expansion rate H- m t during inflation and the new energy scale 
M c . We show how and why different choices of the initial vacuum state (stemming from different 
assumptions about trans-Planckian physics) lead to fractional differences which depend on different 
powers of ao- As we emphasize, the power in general also depends on whether one is calculating the 
power spectrum of density fluctuations or of gravitational waves. 
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I. INTRODUCTION 

The exponential expansion of space in inflationary cos- 
mology leads to the stretching of modes which were in 
the trans-Planckian regime at the beginning of inflation 
into the observable range. This leads to the possibility, 
first raised explicitly in [l| , that trans-Planckian physics 
might be observable today in the cosmic microwave back- 
ground. In earlier work [2, 13 we addressed this issue in 
a simple toy model obtained by replacing the linear dis- 
persion relation of the cosmological fluctuations by new 
dispersion relations which differ from the linear one on 
length scales smaller than the Planck length (the same 
dispersion relations had been used earlier 0, Q in the 
context of an analysis of possible trans-Planckian effects 
on black hole radiation). We were able to construct dis- 
persion relations which give rise to large (order one) cor- 
rections to the usual spectrum of fluctuations, but the 
price to pay is a fine-tuning of the parameters describing 
the model and/or a back-reaction problem. This ques- 
tion has been further analyzed in many papers (see for 
instance Refs. @ SH H Efl It was found that in 

order to obtain measurable differences in the predictions, 
non-adiabatic evolution of the state on trans-Planckian 
scales is required. 

In another line of approach to the trans-Planckian 
challenge to inflationary cosmology, the possibility of 
measurable effects of trans-Planckian physics on ob- 
servables such as CMB anisotropies and power spec- 
tra of scalar and tensor metric fluctuations was studied 
0E0HHE0 in models where the trans- 



Planckian physics is based on stringy space-time uncer- 
tainty relations. In particular, the authors of [l4j found 
a spectrum with oscillations of amplitude ao = H m { /M c , 
where H iii{ is the Hubble parameter during inflation and 
M c a characteristic scale at which the trans-Planckian 
physics shows up, superimposed on the usual scale- 
invariant spectrum, whereas the authors of [l5j found 
only much smaller effects. 

The trans-Planckian problem was also tackled in the 
framework of non-commutative geometry in Ref. [l9l | . It 
was found that the effect is of order (Hi n {/M c ) 4 . It was 
also shown in this article that non-commutative geometry 
implies the presence of a preferred direction which would 
result in a correlation between different multipoles Ce 
and Cg±2- 

In yet another approach to the trans-Planckian issue, 
Danielsson [2(J (see also Ref. |2l|) suggested to replace 
the unknown physics on trans-Planckian scales by assum- 
ing that the modes representing cosmological fluctuations 
are generated mode by mode at the time when the wave- 
length of the mode equals the Planck length, or more 
generally when it equals the length £ c associated with 
the energy scale M c of the new physics which sets the 
initial conditions. There is a one-parameter family of 
vacuum states (a vacua) of a free quantum field in de 
Sitter space which can be considered, and for nontrivial 
a vacua Danielsson found effects of the choice of the ini- 
tial state which are of linear order in the ratio ao 1 , and 
such effects could be seen in observations ^31 2 . Sim- 
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1 Recently, Ref. |22| has shown that effects of the order ao also 
occur in models of power-law inflation. 

2 There has been a concern that nontrivial a vacua are problematic 
from the point of view of interacting quantum field theory l'24l 
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ilar results were found by Easther et al. |28(, whereas 
Niemeyer et al. |2^| have argued that if the modes are 
set off in the initial adiabatic vacuum state when their 
wavelength equals i c , then the effects are of order Cq and 
hence (in usual models) completely negligible. 

U sing an effective field theory method, Kaloper et 
al. [30| have argued that the effects of trans-Planckian 
physics on cosmological fluctuations should be at most 
of the order a 2 , assuming that the modes are in the 
adiabatic vacuum state when the wavelength is equal to 
the Hubble radius (see Ref. [3jj for a criticism of impos- 
ing initial conditions at Hubble radius crossing, and see 
Ref. [13 for counterexamples to the claims of Ref. |30l|'). 

In this paper, we re-consider the calculation of the 
spectrum of cosmological perturbation in the minimal 
trans-Planckian setting | 2C| when mode by mode the ini- 
tial conditions for the mode are set when the wavelength 
equals the Planck length (or, more generally, the length 
scale of the new physics). We find that the overall am- 
plitude of the correction terms (compared to the usual 
spectra) depends sensitively on the precise prescription 
of the initial state, it depends on whether one is studying 
power-law or slow-roll inflation, and it also depends on 
whether one is computing the spectrum of scalar met- 
ric fluctuations or of gravitational waves. Some of the 
"discrepancies" between the results of previous studies 
is due to the fact that different quantities were calcu- 
lated in different models. We show that when the initial 
state is chosen to be the instantaneous Minkowski vac- 
uum, then the deviations of the power spectrum from 
the usual result are of the order erg, in agreement with 
what was found in [2^. In an arbitrary a— vacuum, the 
choice of the value of a has an effect on the amplitude 
of the fluctuation spectrum which is not suppressed by 
any power of Co- However, if a is independent of k, the 
effect will not be distinguishable from a slight change in 
the underlying parameters of the background inflationary 
model. However, in general (and specifically in the choice 
of the vacuum made in 20] , the amplitude of the correc- 
tion term in the power spectrum will have a k-dependent 
(and hence observable) piece which is first order in ao, at 
least in the case of the spectrum of gravitational waves. 

While this paper was being finalized, three preprints 
appeared which investigate related aspects of the trans- 
Planckian problem. In Ref. [3j|, the choice of various 
initial states was related to the minimization of different 
Hamiltonians. In Ref. 34], the predictions of inflation 
for the power spectrum of fluctuations was studied for a 
two parameter class of initial states, and the amplitude of 
the corrections compared to the usual results was seen to 
depend sensitively on which state is chosen. In Ref. [35| , 
the fact that the definitions of the Bunch-Davies and of 
the adiabatic vacua have some intrinsic ambiguities in a 



l25l l2fij . However, very recently it has been shown l27| how to 
define interacting quantum field theory about an a vacuum. 



Universe with a de Sitter phase of finite duration was 
analyzed. 



II. COSMOLOGICAL PERTURBATIONS OF 
QUANTUM-MECHANICAL ORIGIN 

A. General considerations 

The line element for the spatially flat Friedmann- 
Lemaitre-Robertson- Walker (FLRW) background plus 
small perturbations can be written as |36(: 

ds 2 = a 2 (r]){-(l-2(f>)drf + 2(d t B)dx l d-n 

+ [{l-2Tjj)5 t3 +2d l d 3 E + h zj ]dx l dx : >}, (1) 

where the functions cj>, B, ip and E represent the scalar 
sector whereas the traceless and transverse tensor hij 
[hi 1 = hij'i = 0), represents the tensor sector, i.e. 
the gravitational waves. The time r\ is the conformal 
time and is related to the cosmic time t by the relation 
dt = a(rj)di]. It is convenient to introduce the back- 
ground quantity 7(77) defined by 7 = —H/H 2 , where 
a dot means differentiation with respect to cosmic time 
and H is the Hubble rate, H = a/ a. Using the confor- 
mal time we may rewrite 7 as 7 = 1 — Ji' /TC 2 , where 
TL = a 1 / a, and a prime denotes differentiation with re- 
spect to the conformal time. For power-law scale factors, 
a(rf) oc (—r/) q , the function 7(77) turns out to be a con- 
stant, 7 = (g+ 1) /q, which vanishes in the particular case 
of the de Sitter space-time characterized by q = —I. The 
perturbed Einstein equations provide us with the equa- 
tions of motions for the cosmological perturbations. At 
the linear level, each type of perturbations decouple and 
we can treat them separately. 

In the tensor sector (which is automatically gauge- 
invariant) we define the quantity /i T for each mode k 
according to 
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hij (V, x) = \ t^tj £ J dkp& (kK (t?, k)e* kx , (2) 

where pfj (k) is the polarization tensor. The plane waves 
appear in the previous expression because the space-like 
sections are flat. At linear order, gravitational waves do 
not couple to matter and, as a consequence, the equation 
of motion is just given by the per turbed vacuum Einstein 
equation. Explicitly, it reads |37| : 

<+(V-^)// T = 0. (3) 

This equation is similar to a time-independent 
Schrodinger equation with an effective potential given by 
U T (r)) = a" I a. It can also be viewed as the equation of 
a parametric oscillator whose time-dependent frequency 
is given by lo t = k 2 — a" /a, see Ref. [38|. In the vacuum 
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state, the two point correlation function of gravitational 
waves reads 



|^(77, x)/i«fa,x + r)|0> = 



dfc sin fcr 
k kr 



k 3 P h (k). 



(4) 

The power spectrum k 3 Ph(k) which appears in the ex- 
pression for the two point correlation function is given 
by 



k 3 P h (k) 



2fc3 



/'t 



(•5) 



This quantity is a priori time and wavenumber dependent 
but for super- horizon modes, it turns out to be time- 
independent because the growing mode is given by fi T oc 



Let us now turn to scalar metric (density) pertur- 
bations. The two most important differences with the 
gravitational waves are that the scalar sector is gauge- 
dependent and that the scalar perturbations of the met- 
ric are coupled to the perturbations of the stress-energy 
tensor describing the matter. Scalar perturbations of 
the geometry can be characterized by the two gauge- 
invariant Bardeen potentials $ = <f> + (l/a)[(B — E')a]' 
and = tp — H(B — E') [3Jj. From now on, we restrict 
ourselves to the case where the matter is described by 
a scalar field: (p = <fo(ri) + 951(77, x). Fluctuations in 
the scalar field are characterized by the gauge-invariant 
quantity 5ip = ipi + (p' (B — £"). The full set of the scalar 
perturbed (gauge-invariant) Einstein equations is 



J 



dtpo 



di&iQ - ¥) = , (2H' + H 2 )t> + Ht>' + + 2W - -d k d k {<f> - *) = - 



dV 



-(<p' ) 2 $ + cp' 6p' -a 2 ^-5<p 



(6) 
(7) 



where k = Sir/m 2 



being the Planck mass. At this 



point, one has to be careful because the case <p' = 0, 
which corresponds to the de Sitter space-time, is partic- 
ular. In this situation, the solution of the above system 
of equations is simply $ = 0, i.e. there are no density 
perturbations at all (it is also necessary to require that 
the Bardeen potential is finite at infinity). If yd 7^ 0, 
then everything can be reduced to the study of a single 
gauge-invariant variable (the so-called Mukhanov-Sasaki 
variable) defined by 01 



(8) 



It turns out to be more convenient to work with the vari- 
able /i s defined by fx s = —\/2kv. Its equation of motion is 
very similar to that of the gravitational waves and reads 



very 

m 



M s 



(«V7) 



Ms =0 



(9) 



Some remarks are in order here. Firstly, let us stress 
again that the above equation can be used only if ip' ^ 
or 7 7^ because in course of its derivation, we have 
divided by ip' . Secondly, the effective potential for den- 
sity perturbations, U s — (oVtO'V^VtO is different from 
the potential for gravitational waves. Equivalently, its 
time-dependent frequency is also different and given by 
uj s = k 2 — (0^/7)" /(a^/7). Thirdly, in the case of power- 
law inflation, the function 7 is a constant and we have 
U s = U T . Fourthly, density perturbations are often char- 
acterized by the so-called conserved quantity Q defined 



by C = (2/3)(H" 1 $ / + $)/(l + uj) + where u is the 
equation of state parameter, i.e. the ratio of pressure 
to energy density of the background scalar field. This 
quantity is directly linked to fj, s through the equation 
fi s = — 2a % /7C- The integration of Eq. © leads to the 
primordial spectrum of this quantity, namely 



fc 3 P c (fc) 



8tt 2 



Ms 



(10) 



The fact that this quantity is meaningless for the de Sitter 
case, 7 = 0, is obvious. 

Let us now consider a different situation: there are no 
cosmological fluctuations anymore but just a test scalar 
field x( 7 ?j x ) living in a FLRW universe. If we Fourier 
expand the scalar field according to 



1 



1 



dkju(?7, k)e l 



a{r]) (2tt) 3 / 2 

then the Klein-Gordon equation takes the form 



/' 



a 
a 



H = 0. 



(11) 



(12) 



We recognize the equation of motion of gravitational 
waves. The two-point correlation function of the scalar 
field can be written as 



Ix(f7,x)x(?7,x + r)|0) 



dfc sin kr 
k kr 



k 3 P x {k). (13) 



The quantity k 3 P x {k) which appears in the expression for 
the two point correlation function is the power spectrum 
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of the scalar field and is given by 



k 3 P x {k) 



2^2 



(14) 



We also see that the equation of motion for density per- 
turbations is equivalent to the equation for a scalar field 
if the function 7 is a non- vanishing constant. This is 
why, for this class of models, the study of density pertur- 
bations is in fact equivalent to the study of a scalar field. 
However, there is an important exception to this prescrip- 
tion: the de Sitter space-time since, in this case, 7 = 
so that the equation of motion is not given by Eq. @ 
as already mentioned above. Therefore, it is inconsistent 
to first assume that 7 is a constant, then to use Eq. I|1U|) 
in order to calculate the power spectrum and finally to 
particularize the result to the de Sitter case. This proce- 
dure can only lead to the determination of the spectrum 
of gravitational waves and/or of a test scalar field but 
not of density perturbations in the de Sitter space-time. 



B. Power- law inflation 

The case of power-law inflation, for which the scale 
factor can be written as a(rj) — £o(—ri) q , is important 
because the equation of motion for the cosmological per- 
turbations can be solved explicitly. The parameter to is 
a length since we have chosen to work with a dimension- 
full scale factor. Since the function 7(77) reduces to a 
constant, the effective potential for density perturbations 
simplifies to a" /a and becomes identical to the gravita- 
tional waves effective potential. Its explicit form reads 
Z7 S (77) = U T (rj) = q(q — l)/i] 2 . Then, the exact solution 
for the variables ^ s , [A T and \x reads 



Ms 



(k7 1 ) 1 / 2 [A 1 (k)J q _ 1/2 (kr ] ) + A 2 (k)J_ q+1/2 {k V )] . 

(15) 

In the above expression, J v is a Bessel function of order 
v. The two fc-dependent constants A\(k) and A 2 (k) are 
fixed by the initial conditions. 

We now consider the standard calculation of inflation- 
ary cosmology. In this case, the initial conditions are 
fixed in the infinite past. When kr) — > —00, the mode 
function tends to a plane wave with positive and nega- 
tive frequency. The usual procedure requires that 



lim 

k/(aH)^+oc 



Ms 



Ay/n e 



-ikrj 



2k 



(16) 



Let us notice that the origin of the minus sign for 
density perturbations is the minus sign in the relation 
/i s = — \/2kv, v being the quantity which is canonically 
quantized, i.e. which behaves as e~ lkri /V2k in the ultra- 
violet limit. Then, the constants Ai(k) and A 2 (k) are 
completely specified and are given by 



A 2 (k) 



a in(q-l/2) 



(17) 



A 2 {k) = T— , , fc' 1/2 . 
m pl cos(7r<7) 



(18) 



This implies that the mode functions fj, s T can be ex- 
pressed as 

MsT^M) = T^(-^) 1/2 ^ 9/2 <L 9 (-M , (19) 
m pi 

where IfW is the Hankel function of first kind 3 . This re- 
sult allows us to calculate the power spectra. For density 
perturbations and gravitational waves, one respectively 
finds, on super-Hubble scales 



k 3 P c (k) 



t% 7T7 



(21) 
(22) 



where £ rl — W pl is the Planck length and the function 
/(<?) is given by 



/(<?) = - 

7T 



r(i/2-«) 



2" 



(23) 



In the above definition of the function f(q), T denotes 
Euler's integral of the second kind. For the de Sitter 
case, one has f(q = —1) = 1. However, for this case 
the amplitude of density perturbations blows up since 
j(q = —1) = while the amplitude of the gravitational 
waves spectrum remains finite. The origin of the singular 
limit for density perturbations is again the factor at 
the denominator of the expression for the scalar power 
spectrum. In fact, this case must be analyzed separately 
and one can show that density perturbations do not exist 
in de Sitter space-time. We see that power-law inflation 
leads to 

k 3 P h {k) = A T k n t , k 3 P c (k) = A s fc"s- 1 (24) 

with n s — 1 — n T . Explicitly, one has 

n s =2q + 3. (25) 

The case q — — 1 leads to a scale- invariant spectral index, 
namely n s = 1 and n T — 0. 

The power-spectrum of the scalar field can be also cal- 
culated and the result reads (contrary to the previous 
power spectra, the power spectrum of the scalar field is 
a dimension-full quantity) 



k 3 P x (k) 



1 



/(<z)fc 29+2 . 



(26) 



Using the asymptotic behavior of the Hankel function 



11 TTZ 



i(z — 7ri'/2 — n/4) 



(20) 



one can easily check that the mode function fJ. B J , ^ d (k, 77) has 
deed the required limit given by Eq. (16). 
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It is also convenient to expressed this power spectrum in 
terms of the Hubble parameter during inflation H ul { — 
H/a — —q/[£o(—v) 9+1 }- This permits to replace the scale 
£q in the above equation and leads to 



(2n) 2 



(27) 



Of course, in spite of the fact that the time depen- 
dence now appears explicitly, the power spectrum re- 
mains a time independent quantity. In the de Sitter case, 
q = — 1, one recovers the result often cited in the liter- 
ature, namely k 3 P x = [Hi n { / (2tt)] 2 . The spectrum be- 
comes scale-invariant and its amplitude remains finite. 
The situation, except for some unimportant numerical 
factors, is very similar to that of gravitational waves. 

To conclude this subsection, let us emphasize that the 
previous discussion has shown that the analogy between 
density perturbations and a scalar field must be used 
cautiously. The spectral indices are similar even in the 
de Sitter limit but the amplitudes differ radically in this 
limit. 



that the scale factor behaves like a(rf) cx (— r/) _1_e . In- 
terestingly enough, the effective power index at leading 
order depends on e only. At this point, one should make 
the following remark. Using the relation 7 = (q + 1)/<Z, 
the scale factor of power-law inflation can be re-written 
as a(rj) cx (— rj)~ 1 ^ 1 ~~ 1 \ From this expression, one might 
be tempted to write that, in the slow-roll framework, 
0(77) oc (— rj)~~ VC 1-6 ) since 7 = e in this approximation. 
Of course, this expression is inconsistent because it con- 
tains an infinite number of power of e. What should be 
done is to expand the expression —1/(1 — e) in terms of 
e from which we recover that the scale factor is given by 
a{rj) oc (— ■q)~ 1 ~ e . As long as one decides to keep high 
order terms, the whole hierarchy of slow-roll parameters 
should enter the game. This shows that the slow-roll ap- 
proximation does not only consist in naively expanding 
the power index of the scale factor in powers of e. This 
conclusion is reinforced by a study of the cosmological 
perturbations within this approximation. 

The effective potential of density perturbations can be 
calculated exactly in terms of the slow-roll parameters. 
The result is: 



C. Slow-roll inflation 

The slow-roll method is an approximation scheme 
which allows us to go beyond the simple power-law so- 
lutions considered in the previous section. It permits to 
treat a more general class of inflaton potentials. At lead- 
ing order, the approximation is controlled by the slow-roll 
parameters (see e.g. Ref. f° r a new se t of slow- roll 
parameters with a very nice interpretation, see Ref. |43[l 
defined by: 



_Jf_ 

Hip 



2 



V 



H 
H 2 
e-8 



= 1 - 



H 2 



2He 



H 



(28) 



The quantity V(<p) is the inflaton potential. We see that 
7 = e, where 7 is the function that has been introduced 
before. The slow-roll conditions are satisfied if e and 8 
are much smaller than one and if £ = C(e 2 , S 2 , ed). Since 
the equations of motion for e and 8 can be written as: 

±=2e(e-S), A = 2e(e - <5) - £ , (29) 

it is clear that this amounts to considering e and 8 as 
constants if one works at first order in the slow-roll pa- 
rameters. This property turns out to be crucial for the 
calculation of the power spectra of cosmological pertur- 
bations. For power-law inflation the slow-roll parameters 
satisfy: e = 8 < 1, £ = 0. 

The slow-roll approximation can be viewed as a kind 
of expansion around the de Sitter space-time. Indeed, at 
leading order, one has aH w — (1 + e)/r) which implies 



U s fa) = a 2 H 2 [2 - e + (e - 8) (3 - 8) + £] 



(30) 



We have seen before that, in the slow-roll approxima- 
tion, a 2 H 2 ps r]~ 2 (l + 2e). This implies that the effective 
potential reduces to 



*7 s (7j)« A(2 + 6e-3<5). 

T] 2 



(31) 



Since, at leading order, e and 8 must be seen as constants 
in the slow-roll approximation, the equation of motion is 
of the same type as in power-law inflation and the solu- 
tion is expressed in terms of Bessel functions according 
to: 

Ms = (fci7) 1/2 [BiJ-_ 1/a (fc77) + B 2 J_, +1/2 (k V )}. (32) 



The parameter q s appearing in the order of the Bessel 
function is given by 



■1 -2e+S 



(33) 



A comment is in order here: The potential U s depends on 
the scale factor and its derivatives only and we have seen 
before that the scale factor behaves as a(rj) oc (— rj)~ 1 ~ t . 
Therefore, one might think that U s should depend on 
e only. This is not the case. The reason is that U s 
contains terms like e/e (for instance) which are linear 
in 8. First one must calculate all derivatives, replace 
them with their expression in terms of e and 8, and only 
then consider that the slow-roll parameters are constant. 
For gravitational waves, the same lines of reasoning can 
be applied. The effective potential can be written as 
U T (rj) — a 2 H 2 (2 — e) and gives in the slow-roll limit 



2 + 3e 



(34) 
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Therefore, the solution of /x T is similar to the one given in 
Eq. I|32l) , where the effective index of the Bessel function 
is now given by: q T = — 1 — e. This time, the spectral 
index only depends on e as expected from the shape of the 



tensor effective potential. Fixing the initial conditions in 
the infinite past in the standard manner, we arrive at the 
following expressions for the power spectra 



J 



P c (k) = 



V h (k) = 



H 2 



16H 2 



Trm 



2 



1 - 2(C + l)e - 2C(e - 5) - 2(2e - S) In ( — 



1 - 2(C + l)e- 2eln ( — 



(35) 
(36) 



where C = 7 E +ln 2-2 - -0.7296, 7 E - 0.5772 being the 
Euler constant. All quantities are evaluated at Hubble 
radius crossing. The scale fc* is called the pivot scale, 
see Refs. an d E3- One sees again that the de Sitter 
limit e — » is ill-defined for scalar metric fluctuations, 
whereas it is well-defined for gravitational waves. 



that the mode is placed in an a- vacuum according to 



Ms, T (»7fc) = T 



(38) 



/ / \ ,. /w SiT (r? fc )40i : (cfc-4) 
M S , T (%) = ±n/ — ^ , (39) 



III. MINIMAL TRANS-PLANCKIAN PHYSICS 

A. General description 

We now consider the cosmological perturbations of 
quantum-mechanical origin in the framework of the min- 
imal trans-Planckian physics. The new ingredient con- 
sists in assuming that the Fourier modes never penetrate 
the trans-Planckian region. Rather, the main idea is 
that a Fourier mode is "created" when its wavelength 
becomes equal to a new fundamental characteristic scale 
i c = (2ir)/M c . Then, the evolution proceeds as usual 
since the equations of motion of the mode functions n s T 
are taken to be unmodified and still given by Eqs. © 
and @ . In this way, the changes are entirely encoded in 
the initial conditions and there is no need to postulate 
some ad-hoc trans-Planckian physics. The time of "cre- 
ation" of the mode of comoving wavenumber k, rjk, can 
be computed from the condition 



(37) 



which implies that r/k is a function of k. This has to 
be compared with the standard calculations where the 
initial time is taken to be rjk = — oo for any Fourier mode 
k (in a certain sense, this initial time does not depend on 
k). The situation is summarized in Fig.^ 

In the framework described above, a crucial question 
is in which state the Fourier mode is created at the time 
?7fc. There is now no asymptotic region in the infinite past 
where the standard prescriptions can be applied. In this 
article, we consider the most general conditions, namely 



where ui ST = ^Jk 2 — U s T is the effective frequency for 
density perturbations and gravitational waves. The co- 
efficients Ck and dk are a priori two arbitrary complex 
numbers satisfying the condition \ck\ 2 — \dk\ 2 = 1- The 
instantaneous Minkowski state corresponds to Ck = 1 and 
dk = 0. If, in addition, we take r/k — — oo, we recover the 
standard choice, see Eq. since u> 3 T (r/k — — oo) = k. 
We are now in a position where we can compute the mode 
functions and the corresponding power spectra for den- 
sity perturbations and gravitational waves. 



B. Calculation of the mode function 

In this section, we write the scale factor as a{rf) oc 
(— r/) p , where p is a generalized index defined by 

f 



, power-law inflation . 
1-7 ' (40) 



slow-roll inflation 



Since the initial conditions given by Eqs. (|38|l and l|39|l 
are different from the standard ones, we will obtain dif- 
ferent mode functions. The new mode functions that we 
want to calculate can be expanded according to 

ms.tW = as, T (*0/C nd M + PsA k MT d W* > ( 41 ) 

where /ig t ^ nd (?7) denotes the mode functions obtained in 
the standard situation and given by the expression 



2in 



<r d ( fc > ??) = T — (-v) 1/a e-^'*HM_ v (-kr,) . (42) 
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FIG. 1: Sketch of the evolution of the Hubble radius vs time comparing how the initial conditions are fixed in the standard 
procedure and in the framework of the minimal trans-Planckian physics. In the standard procedure, the initial conditions are 
fixed on a surface of constant time 77 = rji with r/i — > —00. In this sense, the initial time does not depend on the wavenumber. 
On the contrary, in the minimal trans-Planckian physics, the initial time depends on the wavenumber and the Fourier modes 
are "created" when their wavelength is equal to a characteristic length. They never penetrate into the trans-Planckian region. 
This causes a modification of the power spectrum which is k dependent since the modification is not the same for all Fourier 
modes as is apparent from the figure. 



In this equation, we have introduced another generalized 
index v defined by 



1 



1 



l- 7 ' 
g s =-l-2e 



k <7t 



-1-e, 



power-law inflation , 
slow-roll inflation , 

(43) 

slow-roll inflation . 



For power-law inflation the generalized index is the same 
for density perturbations and for gravitational waves 



(hence we do not need to distinguish them in the above 
definition). For slow-roll inflation, as already mentioned 
above, they differ. 



The coefficients a ST (fc) and /3 ST (fc) are readily ob- 
tained using the initial conditions given in Eqs. Ij38|l and 



a ST (k) 
f3 s , T ( k ) 



\{c k +d k )e^' 2 ^- 



-2u s (r} k )r)k 



krj h 



H. 



(2) 

Z/2-v 



- H 



(2) 



-2w s _ T (?7fe)?7fc 



krjk 



H. 



(i) 

3/2-v 

r 



- H 



1/2-v 
(i) 

-1/2-v 



2 .Ck — _& w s (r) k )r] k 



l + 2i 



Ck + d k 
Cfe — d k 



H 



(2) 

1/2-v 



Cfe 



w s , x (?7fe)% 



H 



(i) 

1/2-v 



,(44) 



(45) 



where the Hankel functions are evaluated at ~krj k . The alent to the knowledge of the modified mode function, 
knowledge of the coefficients a s T (k) and /3 S T (fc) is equiv- see Eq. I|41|) . The quantity — krj^ can be written as 
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—kr]k = —p/ok where Uk = H(Vk)/M a . The explicit 
expression of 07. can be easily derived. Writing that 
H = p/(arf), we find that <ik — fo(?7fc/?7o)~ p_:L where 
the time 770 is a given, a priori arbitrary, time during in- 
flation which, and this is the important point, does not 
depend on k. The quantity <jq is defined by o~o = Hq/M c . 
Using that i]k = pj (fccfc) and r/o = pj (o~oaoM c ) we finally 
arrive at 



OqM c 



-l-i/p 



(46) 



In the case of slow-roll inflation, the previous expression 
reduces to a> = °~a[k / (aoM c )]~ e . For the de Sitter case, 
p = — 1 and e — 0, the Hubble parameter and therefore 
<7fc are constant. This means that % oc 1/k as expected. 
An important remark is that 07. only depends on p. This 
means that, in the slow-roll approximation, Ok only de- 
pends on the slow-roll parameter e and not on the slow- 
roll parameter S. This is due to the fact that the calcu- 



lation of (Jfc only depends on background quantities. The 
slow-roll parameter S appears in the calculation through 
the generalized index v, i.e. at the perturbed level. We 
now discuss the spectrum obtained according to the ini- 
tial state chosen. 



C. Instantaneous Minkowski vacuum 

In this section, we assume that the initial state is such 
that Cfc = 1 and dk — 0. As a warm up, let us derive 
the spectrum in the particular case of de Sitter, p = 
— 1. As mentioned above, this can only been done for 
gravitational waves since the de Sitter limit is singular 
for density perturbations. For p = —1, we have Ok = o~q 
and —krik = 1/co- The spectrum can be determined 
exactly because the Bessel functions in Eqs. <|44ll and 
(|45|l reduce to ordinary functions in the case p = — 1. 
The result reads 



a T (k) 



: e -Wa-iM)( 1 _2 CT| 



2\-l/4 
07 



i + <r - + (i + <7 Q )tJi — 2g\ 



-i + o-Q + iul + — <ro)yl - 2cto 



The modulus of the two previous expressions can be easily derived from the above relation 



1 (2- 2al-a\ 



1/2 



'0 u o 



1/2 



(47) 
(48) 

(49) 



which coincides with the result of Ref. [2^|. One can check that \a T | 2 — |/3 T | 2 = 1. In addition, one notices that the 
result is independent on k as expected for the de Sitter space-time. Then, the power spectrum on super-horizon scales 
is given by k 3 Ph(k) = (2k 3 / w 2 )\(a T p^ and + /3 T /i^ tand *)/a| 2 where the super-horizon standard mode function can be 
expressed as /x^ and ~ —4iy/Tr(2k)~ 1 / 2 /(m pl kri). The result reads 



k 3 P h (k) = 



V 1 - 2^ 



? 1 a <? ( 2 \ 1 4 (2 
1 - °q - ^0 + CT o sm — + -cr n cos — 



(50) 



r 



This spectrum is represented in Fig.|2 So far, no approx- 
imation has been made. If we now assume that op is a 
small quantity, we can expand the result in terms of this 
quantity. Let us also remark that the terms sin(2/<jo) 
and cos(2/cto) are non-analytic for small values of o~q. At 
leading order, we obtain that |a T | ~ 1 and |/3 T | ~ CTq/2 
in agreement with Ref. (2i|. This leads to the following 
spectrum 



k 3 P h {k) 



1 



Cg sin 



(51) 



Therefore, the corresponding effect in the power spec- 
trum is of order erg, i.e. a tiny effect if Co is small. Note, 



however, that it is not possible to distinguish this effect 
from a change in the parameter H m { of the inflationary 
background. 

We now turn to the case of slow-roll inflation. We 
start with density perturbations. The first step consists 
in calculating the coefficients a s and (3 S . At first order 
in the slow-roll parameters and at leading order in the 
parameter o~q we have 



Jkr/ k 



1 3 r 3 , k 
d em 

2 4 2 a M c 



-ikrj k 



(52) 
(53) 
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0.0 0.2 0.4 0.6 

CT o 

FIG. 2: Amplitude of the gravitational wave power spectrum as a function of the parameter oo in the case where the initial 
state is taken to be the instantaneous Minkowski state. For small values of ao, the leading order correction is cubic in ao- The 
amplitude blows up when cr = 1/V2. 



It is interesting to calculate the modulus of the coefficient 
j3 s . Straightforward calculations lead to 

lAi = Q-3 J -? h ^)'»' (54> 

We see that this quantity depends on both slow-roll pa- 
rameters e and 5, contrary to what was obtained in 
Ref. |2^| where the coefficient (3 S was found not to depend 
on the slow-roll parameter 5. The main reason for this 
discrepancy is that, in Ref. [2j|, the spectrum of a scalar 
field on an unperturbed background cosmological model 
(or, equivalently, the spectrum of gravitational waves) 
was calculated, and not the spectrum of scalar metric 
fluctuations, see the discussion after Eq. lfH)l . If one 
uses the expression of the scale factor a(rj) oc (— r/)~ 1 ~ e 



and inserts this into the expression for the power spec- 
trum for a scalar field on an unperturbed background, 
one does not obtain the correct expression for the power 
spectrum of density perturbations and one misses the 
terms proportional to the slow-roll parameter S. One 
only obtains either the power spectrum of gravitational 
waves or the power spectrum of density perturbations in 
the very particular case of power-law inflation for which 
we have e = S (and for which, in fact, the exact solu- 
tion is known). This is why putting e = S in Eq. I|54|) 
reproduces the result found in Ref. [29|. 

Then, the power spectrum of the conserved quantity £ 
is given by the following expression [to be compared with 



k A Pr 



H 2 I k 

1 - 2(C+ l)e - 2C(e - 5) - 2(2e - S) In — + a%i sin 



Trem 
"3 



An 



:hv 



a Q M c 



j + 3e]n—— 

2 a Q M c 



+ 2(C + l)e + 2C*(e - 5) + 2(2e - S) In 



2 ( k 
— 1 + e + eln— — 

Co V a M c 



-7r(2e — 5) cos 



2 (, k 
— 1 + e + e In 

(T V a Q M 



(55) 
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Thus, the power spectrum has oscillations superimposed 
on the scale- invariant base. Hence, in this case the change 
in the spectrum is "in principle" physically measurable. 
However, the prospects for actually detecting the differ- 
ence of order a 3 between (I55|l and (|35|l are hopeless if the 
parameter ao is small. Indeed if, for instance, we assume 
that do = 10~ 2 , a value consistent with string theory, 
and if we consider that e ~ S ~ 10 _1 then the correc- 



tion to the power spectrum is of order ~ 10~ 7 . There is 
no possibility to detect such a small effect even with the 
Planck satellite. 

For completeness, let us mention the slow-roll result 
for gravitational waves. In this case the coefficients a T 
and (3 T can be deduced from a s and /? s by putting e = 5. 
Then, the power spectrum is given by 



k 3 P, 



i^fl-2(C + l)e-2elnA 





i k Y 




[— (l- 


re + eln — — 




.00 V 


a M c J_ 





-2elnA 



sm 



An 

^0 



:ln- 



a M c 



ne cos 



An 

CO 



:ln- 



a M c 



(56) 



We can also estimate how the consistency check of infla- 
tion is modified. At first order in the slow-roll parameters 
and at leading order in the parameter ao, this quantity 
is the same as in the standard case, namely 

fls ^ = 1& - < 57 > 

This is because at zeroth order in the slow-roll parame- 
ter, the leading contribution in ao is the same for gravi- 
tational waves and density perturbations. 



D. Arbitrary a vacuum state 

We have established that, in the case of the instanta- 
neous Minkowski state, the correction to the power spec- 
tra is of order <7q. It is interesting to see whether this 
conclusion remains true for other initial states. For this 
reason, we now repeat the calculation of the Bogoliubov 
coefficients in the case of a de Sitter space-time with an 
arbitrary a -vacuum state, characterized by a value of c k 
and d k , as the initial state. One finds 



a T (k) = ^(c k +dk)e~™l 2 - l /° a (l-2al)- l l A i + a - ia% + z k {% + a )^l - 2a% 
l3 T (k) = hc k +dk)e^l 2+l '' JO {l-2al)- 1 ' i - i + a + ia 2 Q + z k (i - a ) 0^ 2a 2 



(58) 
(59) 



where we have introduced the quantity z k defined by Zk = {ck — d k )/(c k + dk)- On can easily check that for Ck = 1 
and dk = 0, i.e. z k = 1, the above expressions reduce to Eqs. I)47|l and (|48|) . One can also calculate the modulus of 
these Bogoliubov coefficients, 



\c k + d k \ 2 / ( 1 + z k z* k ) - qg(l + z k z* k ) + crg(l - 2z k z* k ) 
|c fe + d k \ 2 / ( 1 + z k z* k ) -a 2 o{\ + z k z* k ) + - 2z k z* k ) 



+ [z* k + z- i(z* k - z k )a\ 



o] r 



[4 + z k + i(z* k - z k )al} \ , 



(60) 



(61) 



and verify that |a T | 2 — |/3 T | 2 = |c; 12 
spectrum exactly. The result reads 



CfcT — \dk\ = 1- From the two previous relations, one can determine the power 



k 3 P h = 



mH 2 ni \c k + d k \ 2 r 1 



-(l + z k zi)(l-vZ) + -£(l-2z k zi) 



+~(1 - Zir-i, )V-iaZ - 1 ! »«( ) + (J - :/, )a„»m 



Co — Co cos — — 2 sin — 



ao 



ao 



-a*(z* k - z k )Jl - 2a 2 + -a (z* k - z k )(a 2 - 2W1 - 2a 2 cos — 



ao 



(62) 
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As is apparent from the above expression, the correction in the power spectrum compared to the usual results is of 
order unity, i.e. not suppressed by any power of <7q. If one argues that the Bogoliubov coefficients themselves should 
be of linear order in a® then the correction terms would also be linear in <jq. On the other hand, unless the Bogoliubov 
coefficients depend on k, the effect is simply a change in the amplitude of the spectrum, and can hence be absorbed in 
a redefinition of the background. Thus, the effect is not physically measurable. On the other hand, if the coefficients 
depend on k (as discussed e.g. in the analysis of |2lJ |) there is a large measurable effect on the power spectrum. For 
the instantaneous Minkowski state, one check easily that only the two first terms in the curly brackets survives and 
that the corresponding expression reduces to the one found previously. If we expand this spectrum in powers of a a , 
we find that 



k 6 P h 



Ok 



■dA 2 



Zk cos 



CO 



00 



CO 



i(z* k - z fe )cos 



CO 



^0 



(63) 



We conclude that the result obtained in the previous sub- 
section for the instantaneous Minkowske vacuum, namely 
that the corrections to the spectrum are suppressed by 
three powers of Co appears to be very particular to the 
choice of that state. 

However, the previous analysis does not cover all the 
possible cases. Indeed, as we are going to demonstrate, 
the proposal put forward in Ref. [2i| corresponds in fact 
to a case where the complex number Zk can also depend 
on Co- The previous study assumed that Ck and dk were 
pure complex numbers. Therefore, in Danielsson's case 
[2f| . the analysis needs to be redone. 



64tt 



a 2 (r))d 4 x. 



In the first of the two previous expressions, is the 
inverse of the FLRW metric. The next step consists in 
inserting the expression for hij , see Eq. , into the ac- 
tion. This gives 



52 = ~jkf^J ^{^aT^'duh' ) , (66) 



E. Danielsson's a-vacuum state 

We start our analysis with the simplest case, namely 
gravitational waves in a de Sitter background. We treat 
this case in some detail since then these calculations can 
be used to study more complicated situations, for in- 
stance slow-roll inflation. The Einstein- Hilbcrt action 



(64) 



where the fields h s (r],x k ), with s = 1,2, are defined by 
the following expression 



h s ( V ,x k )= * *, /dk/i'faje* 1 **. Mi7i 



afa) (2tt) 3 / 2 



expanded to second order (since we are dealing with first 
order equations of motion) reads 



-g^d^W^ihh)^^ (65) 



16ttG 



Therefore, the action of gravitational waves is equivalent 
to the action of two decoupled scalar fields. Inserting the 
Fourier decomposition of the two scalar fields into the 
expression for the action, one arrives at 



S 2 = 



16tt 



2 

E 

s=l 



R 3 + 



d 3 k «)'*«)'- 



(68) 



We can check that this Lagrangian leads to the correct 
equation of motion: we have (we did not take into ac- 
count the overall constant; here the bar over C means 
that we are considering the Lagrangian in the Fourier 



space) 



5C 

Sfj, 



= -2H«r+2(w 2 -fc 2 x/4)' 



(69) 
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SC 



- = 2GO'*-2w«r 



(70) 



completely fixed and they now depend on the parameter 
ctq. Explicitly, the link can be expressed as 



Therefore, the Euler-Lagrange equations given by the 
well-known expression 



d 
dr) 



SC 



(71) 



reproduce the correct equations of motion for the vari- 
ables (/x*)*. Let us notice that we can also vary the 
Lagrangian with respect to (/x* )* to obtain the equation 
of motion for fi a . 

From the Lagrangian formalism that we have just de- 
scribed, we can now pass to the Hamiltonian formalism. 
The conjugate momentum to /i* is defined by the formula 



SC 



m 

pi 
16tt 



(72) 



In real space, the conjugate momentum II s is defined by 
the expression 



SC 



(73) 



We can check that the two definitions are consistent by 
means of the relation 



(2tt) 3 / 2 



dk<e 4k 



(74) 



Danielsson's boundary condition consists in demanding 
that, at the time of creation r]k, one has the usual relation 
characteristic of the standard vacuum state, namely 



IF( % , x k ) = -ik-^a?{r, k )h s { m ,x k ) . (75) 

Using the Fourier decomposition of the scalar fields h s 
and of their conjugate momenta, this last equation boils 
down to 



G<)'--< = -^m* 



(76) 



stand 



stand 



^stand*y _ ( a '/ a ) At stand* + jfc^st 



stand* 



(77) 



So far no approximation has been made and the previous 
relation is general. We now restrict our study to the de 
Sitter case. Using the explicit form of the mode function 
in this case, one arrives at 



i + 2h] k 



-2ikr/ k 



(78) 



which is exactly the relation found in Ref . .20] . From the 
normalization, we deduce that |a T | 2 = 1 + <7q/4. The 
exact power spectrum can now be determined since we 
explicitly know the coefficients a T and /3 T . Performing 
the standard calculation, one finds 



k'Ph 



■ f 2 

(To sm 



^ (2 
-f-COS 

2 V^o 



This expression should be compared with the correspond- 
ing equation found in the case of the Minkowski instan- 
taneous state, see Eq. (|5U|) . In particular, expanding ev- 
erything in terms of ao, one obtains 



1 — (To sin 



(80) 



i.e., a first order effect instead of a third order effect. 
However, once again note that this effect can be repro- 
duced by redefining the background cosmological param- 
eters, and hence it is not a physically measurable effect. 
The spectrum is represented in Fig. [3J 

Let us now turn to the case of slow-roll inflation. 
The goal is simply to evaluate the ratio (3 T /a T given by 
Eq. H77|) to leading order in the parameter <to and to first 
order in the slow- roll parameter e. One finds 



1% 



2 



1 - e In ■ 



0,qM c 



(81) 



This relation is to be satisfied at the time rj = rjk, and 
implies a link between the coefficients a T and /3 T . Since, 
in addition, |a T | 2 — |/3 T | 2 = 1, the coefficients are in fact 



Of course, for e = 0, one recover the exact result obtained 
previously at leading order in <to- Finally the power spec- 
trum for gravitational waves reads 



k 3 P h 



" 2 



2(C+l)e-2eln 



c 



1 + e + eln — — 

cr V a M c 



l-2((7+l)e 

" 2 



2eln- 



e In ■ 



aoM c 



(To7re cos 
I 



1 + e + e In — — 

(T \ aoM c 



(82) 



Note that since in this case the correction terms lead to spectrum, the effect of the correction terms is in principle 
oscillations in the spectrum about a scale-invariant base 
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FIG. 3: Amplitude of the gravitational wave power spectrum as a function of the parameter ao in the case where the initial 
state is the one singled out by Danielsson's condition. For small values of ao, the leading order correction is linear in ao- 



measurable. Since the amplitude is only suppressed by 
one power of erg, the prospects of being able to detect 
such effects in upcoming experiments are good. 

We now explore the same mechanism but for scalar 
metric (density) perturbations. The action is expressed 
in terms of the variable /i s introduced before. It reads 



reads 



S-2 = \J d 4 x ;// r ^i>,,, 0,,, 
In Fourier space, the action reads 
dry / d 3 kL'u*'- 



("V7)" 2 
V M s 



(83) 



5 2 



MsM E 



The conjugate momentum to the variable fi s is given by 
lip = fj,' and therefore Danielsson's boundary condition 



reads 



-ikfi s 



(85) 



where this relation must be evaluated at the time T] = rjk- 
We note that the term (a'/a)fi s , which was present in 
the case of the scalar field, does not appear for density 
perturbations. This can be traced back to the fact that 
the two quantities that are quantized, /i T and fj, g do not 
have the same action. This means that the link between 
a s and /3 S for density perturbations is not the same as 
for gravitational waves. For scalar metric fluctuations, it 



0s 



(/jStand)/ + ik[l 



stand 



(86) 



"s (A«f and *)' + ife/i§ tand * ' 

This equation should be compared with Eq. J77J). The 

absence of the terms a' /a has important consequences. 
In order to guess what the difference is, we can apply 
the previous equation to the de Sitter case, even if in 
principle this is not allowed (see the discussion above). 
One finds 

Eh = \ e -^k Vk _ ( 87 ^ 

afc i + 2krjk — 2ik 2 rj 2 

Expanding the previous expression in <toi we find that 
the first order terms cancel out and we are left with 
P k /a k ~ -(crg/2)e 2i / CT o. Therefore, the result will be of 
order erg and not of order ao as it was the case for a scalar 
field and gravitational waves. We would have obtained 
a linear correction in the power spectrum if Danielsson's 
condition had been 



Me 



k 

-«-M s 
a 



-ikfx s , 



as it is for gravitational waves, instead of (a' s 
see Eq. JSSJ. 

Let us now evaluate Eq. I|86(l consistently and rigor- 
ously in the slow-roll approximation. One finds 
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As expected, the result is quadratic in <jq. Repeating the for the power spectrum 
standard calculations, one can find the explicit expression 



k 3 P r 



H 2 ( h 
l-2(C+l)e-2C(e-5) -2(2e-5)ln 



-2(2e-<5) In- 2eln 



a M c 



2 



eln- 



l-2(C+l)e-2C(e-S)+e- 



aoM c 



+ crg7r(2e — S) sin 



2 



e In • 



a M c 



(90) 



As in the case of gravitational waves, in principle this is 
a measurable effect since it corresponds to oscillations of 
the power spectrum. Typical power spectra are repre- 
sented in Fig. 01 Since the correction is quadratic in <jq, 
the correction turns out to be extremely small and the 
prospects for the detection of such an effect are not opti- 
mistic even with a high-accuracy experiment like Planck. 

Thus, at least with this choice of the precise form of 
the action |4g, Danielsson's prescription leads to trans- 
Planckian corrections of different strengths for gravita- 
tional waves and for scalar fluctuations. Note that this 
also implies that the consistency relation of inflation will 
be modified. Comparing (|82|l and (|90|l . at leading order 
in Co j one finds 



R = 16e< 1 — (Tq sin 



-5-1 1 



:ln- 



a M c 



(91) 



where the ratio has been evaluated at the pivot scale. As 
expected from the previous considerations, the correction 
is linear in ctq- 



IV. COMPARISON WITH ANALYSES USING 
MODIFIED DISPERSION RELATIONS 

In this section, we consider a case where the trans- 
Planckian physics is described by mean of a modified 
dispersion relation. Although, this phenomenological de- 
scription is different from the one considered above, it 
has already been shown that the corresponding spectrum 
can also possess superimposed oscillations |J 0, • I n 
this section, we calculate the spectrum of a scalar field 
(or gravitational waves) living in a de Sitter space-time 
and we compare this result with the corresponding re- 
sult obtained in the framework described in the previous 
sections, see Eqs. (|50f) and <|79[) . In particular, our goal 
is to determine which parameter controls the magnitude 
of the correction to the standard power spectrum in the 
case of a modified dispersion relation (for the case treated 
before, it was the parameter a$). 

The main shortcomings of describing the trans- 
Planckian physics by a modified dispersion relation is 



that we need to assume something about the physics 
beyond the Planck, contrary to the kind of modifica- 
tion envisaged above. Modifications in the power spec- 
trum are obtained if the WKB approximation is vio- 
lated in the trans-Planckian regime. A typical exam- 
ple where this happens is for the dispersion relation 
uj 2 = k 2 — 2b\\k A + 2&i2fc 6 , which represents 

phys phys phys _ phys 

the first terms of a systematic Taylor expansion. This 
dispersion relation is represented in Fig. [SJ 

The equation of motion for the mode function of a 
scalar field is now given by 



a 



M = 0, 



(92) 



where uj is the comoving frequency given by tu(k,i]) — 
a(?y)u; phyE .[fc/a(?7)]. The difference compared to the previ- 
ous section is that the equation determining the behavior 
of the mode function is modified. The mode functions 
depend on whether w 2 > a" /a or w 2 < a" /a. 

Let us consider an expanding cosmological back- 
ground. There will be four time intervals ("regions"), see 
Fig. El As the Universe expands, the physical wavenum- 
ber will decrease. The evolution starts in Region I, where 
the WKB approximation is valid, the term w 2 dominates 
and the mode function is given by 



1 



VMM) 



-if n ui{k,r)&T 



(93) 



We have chosen the WKB vacuum, i.e. we have retained 
only one branch of the most general solution. At a time 
T]i(k), the mode enters Region II, where the WKB ap- 
proximation is violated, and the mode function can be 
expressed as 



fM n ( V ) = C+(k)a( V )+C^k)a( V ) 



" d7 



a 2 (r) 



(94) 



The coefficients C±(k) can be determined from joining 
the mode function and its derivative at the matching 
point rji(k). After more expansion [at the time 772 (fc)], the 
mode leaves Region II and enters Region III, where the 
WKB approximation is restored. However, now the mode 
function has two branches, the second branch having 
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FIG. 4: Power spectra given by Eq. I|9()^ . The slow-roll parameters are taken to be e = 0.10005 and S — 0.2 which corresponds 
to a tiny red tilt. The pivot scale is chosen to be k, — O.OlMpc" 1 . The time r/o is chosen to be the time at which the pivot 
scale is "created". This corresponds to fe«/ao = M c . The solid line corresponds to ao = while the dotted line is for ao = 0.02 
and the dashed line for ao = 0.03. 



been "generated" in Region II. It is here that the trans- 
Planckian physics modifies the standard result (which 
predicts f3 k = 0). The mode function can be written as 

a k -i /"„ 0((fc,T)dT 

/Lt TTT (77) = — e "t 

VMM) 

The coefficients a k and [3 k can be found in terms of C± (k) 
by matching the solution (and its derivative) at 772 (fc) as 
previously. Finally, in Region IV, the WKB approxima- 



tion is again violated. This region corresponds to the 
usual super-Hubble region. We have 

Miv {r,) = D + (k)a( V ) + D_(k)a( V ) jf . (96) 

The transition between Regions III and IV occurs at the 
time 773 (k), the time of Hubble radius crossing. Our aim 
is to calculate the coefficients D + {k) which determines 
the spectrum of the growing mode. 

Performing the matching at rj^{k), we obtain 



Miv (V) 



a k e 



iU 3 



03 



+ -=^0(77)03 ( a kl3 e- m * + /3 fc73 V f 



" dr 
a*(r)' 



(97) 



'/:f 



where the subscript "3" means that the corresponding quantity is evaluated at rja(k). In the previous equation, we 
have introduced the notation 7^ = u)'/(2ui) + iuo + TL and Q 3 = f v , 3 , w(fc,r)dr. On the other hand, a lengthy but 

straightforward calculation shows that the coefficients a k and 0k are given by 



\/4wiW2 



0-2 * Ol 

—72 7i - aia 2 7i7 2 

ai a 2 



0k = ~ 



\/4wiW2 



a 2 a,i 

— 72 7i 

ai a 2 



7)2 dr 
l a a(r) 
" 2 dr 



aia 2 7i7 2 



-i(n 2 +ni) 



* « 2 (t) 



(98) 
(99) 
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FIG. 5: Typical example of a dispersion relation which breaks the WKB approximation in the trans-Planckian regime but 
allows a non-ambiguous definition of the initial state. 



These expressions are similar to those found in Refs. @ and H3. So far, no approximation has been made. To go 
further, one has to take into account the fact that the difference \i]2(k) — ?7i(fc)| cannot be too large. Otherwise, this 
would mean that particles production in region III is too important and, as a consequence, that the calculation is not 
valid due to this back-reaction problem. Therefore, we can write 771 = 772(1 + A) and perform an expansion of the 
coefficients ak and (3k in terms of the parameter A. The result reads 



Oik 



3 i(n 2 -ni) 



1 + -W2?72 1 



Qi 



A 



-0(A 2 ). 



-i(fi 2 +fii) 



i I ^ 2 



uo 2 a 2 



A + C(A 2 ). (100) 



In this expression, Q is the parameter which controls the accuracy of the WKB approximation. It is defined by 
Q = —lu"/(2uj) + 3(ci/) 2 /(4w 2 ) (3^. As expected, if we send the parameter A to zero, then (3k vanishes and ak 
becomes unity. Inserting the above expressions for ak and (3k into Eq. Ij97(l . one arrives at 



\D + (k)\ 2 = 



1 



2a 2 uj 3 



1 - LU2V2 1 1 — 2 



uj 2 a 2 



'73 



uj(k, r)dr 



in 



A + C(A 2 ) 



(101) 



Notice that in order to establish the previous expression we have just assumed that A is a small number but we 
have not assumed that k c /H or fc*/i? are small. As expected, The correction to |D + (fc)| 2 is given by the Bogoluibov 
coefficient (3k- 

We now need to calculate explicitly the previous quantities. For this purpose, we introduce an approximate piecewise 
form of the dispersion relation considered before, namely 



UJ 



a(k c 



s Khys + (a + l)k c 



(a + s)k t 



kphya < k c , 



(102) 



k , ^ k^ 

phys * 



This piecewise dispersion relation is represented in Fig. |SJ The parameter a controls the slope of the dispersion 
relation in the region where the group velocity does not have the same sign as the phase velocity. The parameter 
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FIG. 6: Approximate form of the dispersion relation considered in Fig. Region II is the region where the WKB approximation 
is not satisfied. 



s controls the slope in the region k phys > fc*. Assuming that the space-time is de Sitter, one can also compute the 
comoving frequency ui = auj vtiys (k/a) 



iv(k,r)) 



a + 1 k 

—ak — 



sk + - 
V 



r) H 
(a + s)— -( a + 1)^2. 



^phys < 1 



^phys > ^* 



(103) 



r 



Obviously, one has to choose k c > afc*/(a+l) in order to 
insure that the frequency remains positive. The various 
times of matching can now be calculated very simply. 
They read 



m(k) 



l 

ks 

1 

k 

1 

k 



V2+(a + s)^-(a + l)^ 



V2 



1 k r 



a 



H 



The times at which > 



= k r , and u> 

rs O p 



(104) 

(105) 
(106) 

= fc», i.e. r/ c (k) 



and ry*(fc) respectively (see Fig. EJ), can also be deter- 
mined easily. They are given by r] c (k) = —{\/k){k c /H) 
and = -{l/k){K/H). 

We have based the calculation of |D+(fc)| 2 on the as- 



sumption that the parameter A is small in order to avoid 
a back-reaction problem. This means that the time spent 
by the modes of interest in the region where the WKB 
approximation is violated is small. In turn, this requires 
a link between the two scales k c and fc* which charac- 
terizes the shape of the dispersion relation and the Hub- 
ble constant H which characterizes the "velocity" with 
which a mode crosses the region where the WKB approx- 
imation is not valid. Using the expressions of r]i(k) and 
772 (&), one finds that the link between k*, k c and A can 
be expressed as 



H 



1 



a a[a + s) 



A 



V2+(a+l)^ 



(107) 



If fc c ^> H and fc* ^ H, then one has fc* ~ (a + l)k c /a, 
as expected. The final result can be expressed in terms 
of k c /H and A only. 
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We can now calculate each term present in Eq. (|101|l . 
Straightforward calculations show that Qjio 2 \i = — (a + 
l)k c /{2V2H) + 3(a + l) 2 k 2 c /(16H 2 ), a"/(aw 2 )| 2 = 1 and 
wry 1 2 = — v2- In order to calculate the integral appearing 
in the argument of the sine function in Eq. (|101|) , one has 
to cut it into several pieces and to use the corresponding 
form of the piecewise dispersion relation. One obtains 



</3 



uj(k, r)dr 



ni 



/ uj(k, r)dr - V2A 

Jr/2 



(108) 



= (V2- 
-V2A. 



(a + 1) - V2H/k c 



(109) 



We have to insert this last expression into the sine func- 
tion and expand the resulting expression once more in 
terms of the parameter A. The final result reads 



fc 3 



£ 1 + V2 



H 

+0(A 2 ) 



1 kr- 



3(a + l) 



2 K 



2^2 H 



16 H 2 



(a+l)^ln 



1 - V2H/k c 



A 



(110) 



As expected the result does not depend on k since the 
scalar field lives in de Sitter space-time. It is worth notic- 
ing that we have not assumed anything about the ratio 
k c /H in order to obtain the previous expression (only the 
parameter A was supposed to be small). The previous 
expression is the main result of this section. It should be 
compared with Eqs. (|50|l and l|79|) obtained previously. 
We see that the magnitude of the correction is no longer 
controled by any power of the ratio of two scales as it was 
before. The magnitude of the effect is determined by the 
time spent in the region where the WKB approximation 
is violated which in turn depends on the shape of the 
dispersion relation that was assumed. 

Finally, let us return to the back-reaction problem. 
There is no back-reaction problem if |/3fc| <C 1. Using 
the expression of (3k derived before and assuming that 
k c /H ^> 1, we see that the parameter A must satisfy 
A <C H 2 /k 2 . This amounts to a severe fine-tuning of 
the scales k c , and of the Hubble parameter H during 
inflation in order to satisfy this condition. 

V. CONCLUSIONS 

We have analyzed the magnitude of correction terms 
to the power spectra of scalar metric (density) fluctua- 
tions and gravitational waves in inflationary cosmology 
under the assumption that fluctuation modes are gener- 
ated when their physical length scale equals some critical 
length determined by the unknown Planck-scale physics, 
but without modifying the equations of motion for the 
fluctuations. The magnitude of the correction terms can 
then be expressed as a function of the dimensionless ra- 
tio (To = Ho/M c , where Hq is the characteristic Hubble 
expansion rate during inflation, and M c is the mass scale 
at which the new physics sets in. 



It is important to realize that the magnitude of the 
correction terms is in general different for gravitational 
waves and for scalar metric fluctuations - a point not real- 
ized in some papers on the "trans-Planckian problem" of 
inflationary cosmology. In addition, the magnitude of the 
correction terms depends sensitively on the initial state 
chosen. We have shown that for the local Minkowski vac- 
uum state, the correction terms are of the order ctq (in 
agreement with the results of pflj). whereas for nontriv- 
ial a- vacua the effects are much larger. If the Bogoliubov 
coefficients which describe the mode mixing do not de- 
pend on Co, then the correction terms can be of order 
unity. In the case of Danielsson's a-vacuum for which 
the Bogoliubov coefficients depend on Co , the corrections 
to the gravitational wave spectrum are suppressed by one 
power of er (in agreement with the results of [2D]), the 
corrections to the scalar metric fluctuation spectrum by 
two powers of cto- Thus, the consistency relation for fluc- 
tuations in inflationary cosmology obtains corrections of 
linear order in ao, as already emphasized in [Tfil ]. 

In the final section of the paper we compared the re- 
sults obtained in earlier sections with the results obtained 
by assuming that trans-Planckian physics leads to a mod- 
ified dispersion relation. In this case, corrections to the 
usual power spectra of fluctuations can be obtained which 
are not suppressed by any small dimensionless combina- 
tions of energy scales in the problem. However, demand- 
ing that the back-reaction remains under control leads to 
severe fine-tuning requirements on such models. 
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